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The mutual dipole-dipole interaction of atoms in a trap can affect their fluorescence. Extremely 
large effects were reported for double jumps between different intensity periods in experiments 
with two and three Ba + ions for distances in the range of about ten wave lengths of the strong 
transition while no effects were observed for Hg + at 15 wave lengths. In this theoretical paper 
we study this question for configurations with three and four levels which model those of Hg + 
and Ba + , respectively. For two systems in the Hg + configuration we find cooperative effects of up 
to 30% for distances around one or two wave lengths, about 5% around ten wave lengths, and, 
for larger distances in agreement with experiments, practically none. This is similar for two V 
systems. However, for two four-level configurations, which model two Ba + ions, cooperative effects 
are practically absent, and this latter result is at odds with the experimental findings for Ba + . 

PACS numbers: 42.50.Ct, 42.50.Ar, 42.50.Fx 



I. INTRODUCTION 



The dipole-dipole interaction is ubiquitous in physics 
and for example responsible for the ever present van der 
Waals force. It is also important for envisaged quan- 
tum computers based on atoms or ions in traps. Con- 
siderable interest in the literature has also been aroused 
by its cooperative effects on the radiative behavior of 
atoms pj. In an as yet unexplained experiment 0, Q 
with two and three Ba + ions, which exhibit macroscopic 
light and dark periods, a large fraction of double and 
triple jumps was reported, i.e. jumps by two or three 
intensity steps within a short resolution time. This frac- 
tion was orders of magnitudes larger than for indepen- 
dent ions. The quantitative explanation of such a large 
cooperative effect for distances of the order of ten wave 
lengths of the strong transition has been found difficult 

IB lil EH Experiments with other ions showed no 
observable cooperative effects 0, , in particular none 
were seen for Hg + for a distance of about 15 wave lengths 
|13| . More recently, an unexpected high number of simul- 
taneous quantum jum ps in a linear chain of trapped Ca + 
ions was reported H3 w hile no such effects were found in 
another experiment |l5j using the same ion species and 
a similar setup. 

Systems with macroscopic light and dark periods 
can provide a sensitive test for cooperative effects of 
the dipole-dipole interaction. These periods can occur 
for multi-level systems where the electron is essentially 
shelved for some time in a meta-stable state without pho- 
ton emission [l6| . For two V systems with macroscopic 
light and dark periods the effect of the dipole-dipole inter- 
action was investigated numerically in Ref. |17| and ana- 
lytically in Ref. |2| and shown to be up to 30% in the dou- 
ble jump rate compared to independent systems. Moni- 
toring the dipole-dipole interaction of two V systems via 
quantum jumps of individual atoms was inves tigate d in 
Ref. 0| . The experimental systems of Refs. |3Tll3l IT^| 
are, however, not in the V configuration so that the re- 



sults of Ref. do not directly apply. 

The experiment of Ref. ^| used two Hg + ions, with 
the relevant three levels as in Fig. which we call a 
D configuration. From a theoretical point of view two 
such systems were studied in Ref. |20j for the special 
case Ai, A3 <C r <C A2, where r is the distance and the 
wave lengths refer to the respective transitions of Fig. 
and for this case no cooperative effects were found. The 
general case will be treated explicitly further below. The 
cooperative effects found here are of similar magnitude 
as for V systems, and for distances of the above range the 
result of Ref. [2(j is confirmed. Our results are also in 
agreement with the experimentally observed absence of 
cooperative effects for distances of about 15 wave lengths 

m 

The levels of Ba + used in the experiment of Refs. 0,0] 
are depicted in Fig. 0(a). The ground state 6 2 Si/ 2 and 
the two upper states Q 2 Pi/2 and 5 2 D 5 / 2 constitute a 
strongly driven fluorescing A system which provides the 
light periods. Only when the system is in the ground 
state can the weak incoherent driving of the 6 2 Si/2 - 
6 2 P3/2 transition populate the meta-stable 5 2 D 5 / 2 state, 
with ensuing dark period. Therefore the details of the 
two upper states of the A system play no significant role 
for the transition to a dark period, and therefore these 
two states are here replaced by an effective single level. 
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FIG. 1: Three-level system in D configuration with fast tran- 
sitions (solid lines) and slow transitions (dashed lines) . 
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FIG. 2: (a) Relevant level scheme of Ba+ HQ. For circled 
levels see text, (b) Effective four-level system for Ba + . Strong 
coherent driving of the |1) — j3) transition by a laser, weak 
incoherent driving of the jl) — j4) transition by a lamp, weak 
decay of level 1 2} . 



This leads to the four-level configuration of Fig. [21 (b). 
The present paper is, to our knowledge, the first to the- 
oretically investigate possible cooperative effects for two 
such four-level systems. Surprisingly, these effects turn 
out to be much smaller than for two V systems for dis- 
tances r > A3, and this shows that cooperative effects 
sensitively depend on how the meta-stable level is pop- 
ulated. Our results for two four-level configurations are 
at odds with the experimental findings of Ref. [3 on the 
magnitude of double jump rates for two Ba + ions pf|. 

The methods presented in this paper can be carried 
over to describe the Ca + experiments of Refs. [Til llj| 
although this would of course require the use of a different 
level system. 

The plan of the paper is as follows. In Section ITU we 
treat two dipole-dipole interacting D systems and explic- 
itly calculate the transition rates between the various 
light and dark periods as well as the double jump rate. 
This is done by means of Bloch equations. In Section 
IIIII the method is carried over to two four-level systems 
of Fig. |21 (b) and the transition rates are calculated. In 
the Appendix a direct quantum jump approach [22| is 
outlined for two D systems. 



II. TWO DIPOLE-INTERACTING D SYSTEMS 

The D configuration, as displayed in Fig. 2] is a model 
of the level system of Hg + used in the experiments of 
Refs. 0,|23j- The transition — 13) is driven by a strong 
laser. Level |3) can also decay via a slow transition to 
the meta-stable level |2). For simplicity all transitions 
are treated as dipole transitions. 

In the following we will investigate two dipole-inter- 
acting D systems a fixed distance r apart and calculate 
the transition rates between the three types of fluores- 
cence periods. This will be done in this section by means 
of Bloch equations. In the Appendix the efficient quan- 
tum jump approach will be applied to two D systems, not 
only confirming the Bloch equation result but also giv- 
ing higher-order terms. For simplicity the laser direction 
will be taken as perpendicular to the line joining the two 
systems. Rabi frequency and Einstein coefficients satisfy 



A 3 ,Q 3 > A U A 2 . (1) 

It is convenient to use a Dicke basis, 

\g) = |1)|1), |e 2 ) = |2)|2), |e 3 ) = |3)|3) 
M = (|*>|i) + |i>|i»/V2, (2) 

K) = (K)b)-|j)N»/V2i. 

In Fig. (a) the level configuration for two D systems is 
displayed in this basis, with the slow decays neglected, 
while in Fig. [21(b) only the slow decays are shown. From 
Fig- El ( a ) it is seen that without slow decays (i.e. for 
A\ = A2 = 0) the configuration decouples into three 
independent subspaces, denoted by So,S\,S2-, with 



50 = {|e 2 )} 

51 = {|S12>, |ai2>, |S23>, |fl23>} 

52 = {\g), |si 3 >, |ai 3 ), |e 3 )} • 



(3) 



By means of the conditional Hamiltonian, H cont x, and 
the reset operation, 1Z, the Bloch equations can be writ- 
ten in the form [24J 
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K(p) 



(4) 



For two D systems one finds by the same method as in 
Refs. II El 
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FIG. 3: Level configuration of two D systems in the Dicke 
basis, (a) Slow transitions omitted, (b) Transitions with rate 
A2 ± Re C2 (dotted arrows) and transitions with rate Ai ± 
ReCi (dashed arrows). Line shifts due to detuning and to 
Im d axe omitted. 
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where A3 is the detuning of the laser. The complex dipole 
coupling constants Cj depend in an oscillatory way on the 
distance r, 
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with fej = 2Tr/\j and #j the angle between the corre- 
sponding dipole moment and the line connecting the sys- 
tems. For maximal effect we take 8j = ir/2 in the follow- 
ing. The real part of Cj leads to changes of the decay 
constants and the imaginary part to a level shift in the 
Dicke basis, as seen from the expression for H conl ±- The 



reset operation can be written in the form 
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To determine the transition rates we write Eq. Q in a 
Liouvillean form as 

p = £p = {£ (A 3 ,C 3 ,n 3 ,A 3 ) + L 1 {A 1 ,A 2 ,C 1 ,C 2 )}p 

(8) 

where the super-operator Ci(A±, A 2 , C\, C 2 ) is a pertur- 
bation depending on the small parameters, and employ 
the following important property of the time develop- 
ment. Starting with an initial state in one of the sub- 
spaces Si, the system will rapidly — on a time scale of Q^ 1 
and A^ 1 — approach one of the quasi-stationary states 
Pssj. Thereafter — for times much larger than f^ 1 and 
A^ 1 , but much smaller than A^ 1 and A^ 1 — small popu- 
lations in the other subspaces will build up until eventu- 
ally, on a time scale of A^ 1 and A^ 1 , the true stationary 
state is approached. Hence we consider a time At with 



A^ 1 ,^ 1 < At < Ay 1 , A^ 1 



(9) 



and calculate p(to + At) for initial p(<o) — Pss,i- 
The quasi-invariant states are easily calculated from 
£ p ss ,i = as 



p ss ,o = |e 2 >(e 2 ], 
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where 



{N-nl(2Al + 3nl + 8A 2 3 )}\g)(g\ 

+ nl(2Ai + n| + 8A|)|si 3 )(a 1 a| 
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As in Ref. [2j one has in perturbation theory 



{A 2 3 + ^A 2 3 ){\C 3 \ 2 



p(to + At) = 



(13) 



At 



p ss ,r+ J dTe C0t C lPs ^+O({A 1 ,A 2 ,C 1 ,C 2 ) 2 ), 



but, unlike Ref. 2], Cip ss ,i is not a superposition of just 
the eigenstates for nonzero eigenvalues of Co but also of 
the Pass's. We therefore decompose £ip ss ,i into a super- 
position of all eigenstates (matrices) of Co, 



£lPss,i = \] 0!ijP Bat j + p 



3=0 



(14) 



where p contains the contributions from the eigenstates 
for nonzero eigenvalues of Co- For later use we note that 
these eigenvalues are of the order of A 3 and il 3 . The co- 
efficients ay can easily be determined by means of the re- 
ciprocal (or dual) eigenstates where only those for eigen- 
value of Co are needed. They are denoted by p* s and 
are defined through 



i,j = 0,1,2 



Ti(pilC A) = for any matrix A. 



(15) 



(16) 



The latter means Clpl s = 0, with the adjoint C\ defined 
with respect to a scalar product given by Tr(A^B). Then 
one has 



ay = Tr(^|£ip SS)i 



The reciprocals p z ss are easily determined as follows. Since 
the Bloch equations conserve the trace one has 



= Trp = Tr£oP 

for any p. Thus 

= Tr(l£ p) = Tr((4l)p) 



(18) 
(19) 



for any p and therefore £gl = 0. Now 1 can be written 
as a sum of terms purely from Sq,S\, and 52 and, since 
the subspaces are invariant under Co, these terms must 
be annihilated by C\ individually. This yields 

Pi = Ie 2 )(e 2 |, 

pl s = |S12)(S1 2 | + \ax 2 )(a 1 2\ + |S23)(S23| + |023)(023|, 

P 2 ss = ItfXffl + |si3>(si 3 ] + |oi3>(ai 3 | + |e 3 )(e 3 | 

since the sum of the right-hand sides indeed yields 1 and 
the normalization condition of Eq. (|15fl is fulfilled. From 
Eq. (|17f) one then obtains the ctij. Inserting now C\p SSt i 
into Eq. 1)1 -ill one obtains 

p(t + At) = 

y / 2 \ 

p ss ,i + / dr 2J aijp ss ,j + e CoT p (20a) 
o V'=° J 

2 

= Psa,i + ^oiijPssjAt + (e - Co)~ x p (20b) 

where for the p term the upper integration limit can be 
extended to infinity since p belongs to nonzero eigenval- 
ues of Co and is therefore rapidly damped. Now, Cq 1 
is of the order of A^ 1 and fij 1 on p, and thus the last 
term in Eq. i|20bjl is of the order of C\/ (A 3 ,£l 3 ) which 
is much smaller than ctij At ~ C\At, by Eq. ©. There- 
fore the last term in Eq. I|20b|l can be neglected, and this 
equation then reveals that the a^ 's have the meaning of 
transition rates from subspace Si to Sj, i.e., 



(21) 



The transition rates are now obtained from Eqs. I|21|l . 
(H3, and (0-011) as 
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One sees that, for two D systems, only p 2 i depends on 
the dipole coupling constants C 3 in first order. This con- 
trasts with two V systems where both p 2 i and pi 2 de- 
pend on C 3 . The physical reason for this is that transi- 
tions between bright periods for two D systems are due 
to decays and not due to the laser. Hence the transition 
rates should essentially be governed by the Einstein coef- 
ficients of these decays on the one hand and by the pop- 
ulation of the states in the initial subsystem on the other 
hand. Therefore the parameter C 3 enters only through 
the quasi-stationary state /? ss ,i of the initial subsystem. 
The absence of a linear C\ and C 2 dependence can be 
understood from Fig. [3] (b) as follows. For most slow 
transitions between two subspaces there are two chan- 
nels with rates Aj ± Re Cj so that Re Cj cancels. States 
with a single decay channel lie in Si and, by symmetry, 
they appear in pairs with different sign of Re Cj . 

From the transition rates pij one obtains the double 
jump rate, tidj, by the formula 2] 

Pl0P0lPl2P21 Arr ro „s 

TlBJ = 2 ; ; AJ D J (23) 

P01P21 + P21P10 + P01P12 

where ATdj is the defining small time interval for a dou- 
ble jump. Significant cooperative effects occur only as 
long as f2 3 and A 3 are at least an order of magnitude 
smaller than ^4 3 . When compared to non interacting 
systems, the cooperative effects are up to 30% for dis- 
tances between one and two wave lengths and 5% around 
ten wave wave lengths, similar as for two V systems. 
For longer distances they are practically absent and this 
is consistent with the experimental results of Ref. [l3j . 
Fig. 0] shows P21 and «dj versus the relative distance 
r/A 3 for typical parameters, where A 3 is the wave length 
of the strong transition. 

Our explicit results for arbitrary r confirm the large- 
distance result of Ref. [2(j who argued that for Ai, A 3 <C 
r <C A2 cooperative effects are "suppressed by the rapid 
decay on the fast transition". In fact, we find in the 
Appendix that these effects are only to first order inde- 
pendent of the coupling parameter C 2 ; the second order 
contributions in C2 are, however, negligible for the ex- 
perimental values of Ref. 01 • 



III. TWO DIPOLE-INTERACTING 
FOUR-LEVEL SYSTEMS AS A MODEL FOR 
TWO Ba+ IONS 

The experiments of Refs. used Ba + . As explained 
in the introduction we model the relevant level scheme by 




FIG. 4: Transition rate P21 and double jump rate mdj for two 
dipole-interacting D systems. The dashed lines show the case 
of independent systems. Parameter values are A\ — Is -1 , 
A 2 = 1 s"\ A 3 = 4 • 10 8 s _1 , Q 3 = 5 • 10 7 s _1 , and A 3 = 0. 



the effective four- level configuration in Fig. [2] (b) . The 
|1) — |4) transition is driven weakly and incoherently by 
a lamp, while the |1) — |3) transition is driven coherently 
by a stron g la ser. This time the Bloch equation can be 
written as |26| 
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{4 {0) (A 2 , A 3 , Ai, Q 3 , A 3 , C 3 ) + 4 (1) (C 2 , C 4 )} P 



where TZw(p) describes the incoherent driving as in 
Ref. j2(J and is given explicitly below. The Dicke states 
are defined in analogy to Eq. (OL and i/cond and 1Z(p) 
can be calculated as in Refs. pi l25l] as 
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The lamp term is obtained as in Ref. (2(| as 

n w { P ) = (27) 

W(R 4+P R\ + + R 4 -pR\_ + R\ +P R 4 + + R\_pR 4 -) 



where W is the product of the spectral energy density of 
the lamp and the Einstein B coefficient of the |1) — |4) 
transition. 

Now the procedure is similar as for the D system. The 
Liouvillean Cq possesses three (quasi-) stationary states 
Pss,o 5 /°ss,i, and p SSi2 which coincide with those for the D 
systems in Eqs. <|10I12[1 and which are associated with the 
dark and the two bright periods. As before, one calcu- 
lates p(to + A£) as in Eq. (| 1 31) and decomposes C 4 p SS: i as 
in Eq. I|14[l . Now, however, the reciprocals p l ss are more 
difficult to determine since |4) can decay into |1) as well 
as |2). An exact solution of Clp l ss = is rather elaborate. 
We therefore decompose 



4 = 4 (0) (A 2 , A 3 , A 4 , 3j A 3 , C 3 ) +4 (1) (C 2) C 4 ). (28) 



and, by Maple, have calculated p\ s to first order in pertur- 
bation theory with respect to C 2 and C 4 , with the same 
constraint as in Eq. i|15|) . The lengthy result will not be 
given here explicitly. The transition rates are again given 
by 



Pij = Tr(p{t£ip ss ,i) 



(29) 



and one obtains for two dipole-interacting four-level sys- 
tems of Fig. [3 (b) to first order in C 2 and C 4 



Poi 
Pio 
P12 



2A 1 



A 2 W{A 2 3 + fi§ 



4A|) 



(A 2 + A 4 )[Al + 2tt§ + 4A§] 

A! 



(30a) 
(30b) 
(30c) 



and 



7 



P21 



2A 2 W^ As 



•4A|)(AI 



2fi§ 



4A|) 



(A§ + 4A|)(|C 3 | 2 + 2A 3 ReC 3 + 4A 3 ImC 3 ) 



( J 4 2 + A 4 )(^ + 2^ 



= 2A 2 W 



4A§ 



{A 2 +A 4 )[A 2 3 + 2n 2 3 + 4Al} 



4A|) 2 + (4| - 
2ReC 3 



-4A§)(|C 3 | 2 + 
A s Ql(Al 



2A 3 ReC 3 
+ 4A|) 



4A 3 ImC 3 



(A 2 + A 4 )[Aj 



4ImC 3 



A3fii(A§ 



2ft 2 . 

4A|) 



•4A§]3 



(A 2 +A 4 )[A2 + 2r!2 + 4A2]3 



0(C 2 ). 



(30d) 



It is seen that poi, Pio, and P12 are independent of the 
coupling parameters and are thus the same as for non 
interacting systems. 

These results for two four-level systems show great sim- 
ilarity with those for the two D systems of the proceed- 
ing section. In both cases only p 2 ± depends to first or- 
der on C 3 , the coupling parameter associated with the 
laser-driven transition. However, cooperative effects are 
significantly smaller for the two four-level systems. For 
fixed laser detuning, the effect of C 3 becomes maximal 
for fi 3 = ly/Al + 4A 2 . For this value of ft 3 , 

Fig. shows the transition rate p 2 i from a double in- 
tensity period to a unit-intensity period and the double 
jump rate tidj over the relative distance r/A 3 , with the 
other parameters as in the experiment [jj 0] . Despite the 
optimal choice of the Rabi frequency, f2 3 , the deviations 
from the value for non-interacting systems are very small. 
Already for a distance of about a wave length A 3 , they 
are not more than 1% for p 2 i when compared to non in- 
teracting systems, while for n DJ they are less than 1% . 



IV. CONCLUSIONS 




We have investigated the effect of the dipole-dipole in- 
teraction for two fluorescing systems with macroscopic 
light and dark periods, first for three-level D configura- 
tions and then for four-level systems. The three-level D 
configuration models the relevant levels of Hg + used in 
the experiments of Ref. and the four- level config- 

uration is an effective model for Ba + , used in the exper- 
iments of Ref. 0, 0| . For these systems one has macro- 
scopic light and dark periods, and their statistics can be 
a sensitive test of the dipole-dipole interaction. We have 
explicitly calculated the transition rates between the dif- 
ferent light and dark periods by employing Bloch equa- 
tions as well as a direct quantum jump approach. From 
the transition rates the double jump rates are obtained. 

For two D systems the effect of the dipole-dipole in- 
teraction is of similar magnitude as for two V system 
investigated earlier and shown to be up to 30 % for 
distances of the order of a wave length of the strong tran- 
sition and about 5% around ten wave lengths, when com- 
pared to independent systems. For longer distances they 



FIG. 5: Transition rate P21 and double-jump rate tidj 
for dipole-interacting four-level systems, with optimal SI3 = 

\ \J \/5 — 1 \J A\ + 4 A| and all other parameters as in the ex- 
periment The dashed lines show the case of independent 
systems. 



are practically absent and this is in agreement with the 
experimental results of Ref. [13| . We have also recovered 
the special case of Ref. [2(j where distances satisfying 
Ai, A 3 <C r <C X 2 were considered and where an argu- 
ment for the non dependence on the dipole coupling con- 
stant C 2 was given. Here we have shown that this holds 
to first order and that the explicitly determined second 
order terms are negligibly small. 

For the effective model of two Ba + systems our re- 
sults yield very small and hardly observable cooperative 
effects for the double jump rate. This is at odds with 
experimental result in Ref. Our method also applies 
to three Ba + ions, but this is more tedious and requires 
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another paper. Also a theoretical investigation of the ex- 
periments with Ca + 0,^1 is possible with our method. 
For this the calculations have to be carried over to a level 
scheme modeling that of Ca + . 

A further conclusion of our work is the observation that 
the magnitude of cooperative effects due to the dipole- 
dipole interaction sensitively depends on how the meta- 
stable level is populated. 

APPENDIX: QUANTUM JUMP APPROACH 
FOR TWO D SYSTEMS 

The procedure will first be explained for a single D 
system which has just two types of periods, light and 
dark ones. From its level configuration in Fig. ^ it is 
evident that the onset of a dark period is preceded by a 
photon from the |3) — |2) transition, with frequency u>2- 
Hence, starting at to = in |1), the probability density 
for the next photon to occur at time t and to come from 
the 1 3) — 1 2) transition is 

wi U2 (t) = A 2 \(3\ e - iH °°^/ n \l)\ 2 (A.l) 

since iJ con H_gives the time development between photon 
emissions |22j. Then its time integral, 

oo 

P^ 2 = J dtw lLU2 (t), (A.2) 
o 

is the probability for the next emitted photon to come 
from the |3) — |2) transition. Now, let the photon rate in 
a light period be denoted by II ■ Then, after each photon 
of the light period the system is reset to the ground state 



Thus pi2 is given by 



In a similar way one obtains pio and p2i- The transi- 
tion rate poi can be directly read off from the no-photon 
probability e~ 2Alt . One obtains the same results as in 
Section ITT1 when one expands in the small parameters. In 
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